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A NUMERICAL METHOD FOR STUDYING THE TRANSIENT BLADE MOTIONS 
OF A ROTOR WITH FLAPPING AND LEAD-LAG DEGREES OF FREEDOM 


By Julian L. Jenkins, Jr. 

Langley Research Center 

SUMMARY 

The equations of motion for the fully articulated rotor system (i.e., rigid blade with 
flapping and lead-lag degrees of freedom) are derived and the procedure used to obtain a 
numerical solution of the nonlinear coupled equations of motion is outlined. The equations 
can be applied to a rotor system with blades of any airfoil section, mass distribution, 
twist, taper, and root cutout. The hinge geometry is variable; however, the restraints 
considered require that the lag hinge be coincident with or outboard of the flapping hinge. 
Two-dimensional stall and compressibility effects can be accounted for. 

Stability boundaries are defined in terms of an absolute limit (i.e., hinge motion 
amplitudes greater than 90°) for varying blade mass constants as a function of the ini- 
tial azimuth position for introducing a disturbance. In addition, a boundary is presented 
which is based on the criterion of neutral spring forces in the forward quadrants of the 
rotor when the blade is released with an initial flapping- angle displacement at an initial 
azimuth angle of 90°. This boundary is thus a measure of the tip- speed ratio at which 
the nonlinear spring forces in the forward quadrants are approximately zero. A com- 
parison of this boundary with one based on a maximum flapping amplitude limit of 15° 
for a 30-ft/sec (9.14-m/sec) vertical gust indicates good correlation. 

INTRODUCTION 

The problem of defining rotor-blade- motion stability and predicting response to 
disturbances has been a subject of study for many years. These studies have produced 
numerous and ofttimes contradictory results, depending on the simplifying assumptions 
made in order to handle the complicated rotor dynamic system. Such assumptions gen- 
erally leave the applicability of the results subject to question, particularly at extreme 
operating conditions. 

Although most investigations have shown the rotor flapwise motion to be very stable 
below tip-speed ratios of 1.0 (see, for example, refs. 1 to 3), there have been indications 
that inclusion of the nonlinear coupling of flapwise and chordwise blade motions might 


introduce marginal stability characteristics of the rotor (see ref. 4). In any event, con- 
tinuing efforts to increase the operating speed of both the pure helicopter and the com- 
pound configurations require a more exacting treatment of the nonlinear dynamics of the 
rotor system inasmuch as it becomes progressively more difficult to prevent excessive 
blade motions at high tip- speed ratios and under maneuvering conditions at moderate 
tip-speed ratios. Also, the transient response of the rotor to gusts is of more importance 
at higher tip- speed ratios because of increased rotor response sensitivity. 

In an effort to obtain further insight into the problems of blade-motion stability 
and to evaluate the possibility of unstable or excessive blade motions as a result of non- 
linear coupling of flapwise and chordwise motions, the rigid-body equations of motion 
for the fully articulated rotor system (i.e., rigid blades with flapping and lead-lag 
degrees of freedom) have been derived and programed for obtaining a numerical solution 
on a digital computer. This report presents a description of the method of derivation of 
the equations and the procedures used to obtain solutions. In addition, the results of a 
study of the blade-motion stability and response characteristics of the fully articulated 
rotor system are presented and discussed. 

SYMBOLS 

a speed of sound, ft/sec (m/sec) 

ap resultant blade- element acceleration about flapping hinge, ft/sec2 (m/sec2) 

a-L resultant blade-element acceleration about lag hinge, ft/sec2 (m/sec^) 

Aq collective pitch angle at blade root, average value of instantaneous blade- root 

pitch angle around azimuth, deg 

Ai,Bi coefficients of -cos i// and -sin respectively, in egression for 9 ; 

therefore, lateral and longitudinal cyclic-pitch angles, respectively, deg 

b number of blades 

B tip-loss factor (assumed equal to 0.97 herein); blade elements outboard of 

radius BR are assumed to have profile drag but no lift 

c blade section chord, ft (m) 

c d o section profile-drag coefficient 
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c l 

C d ,F 

c d,L 

C M D 


c M t 


dD 

e l 

e 2 


e t 


g 

!e 


Ih 


K SF 


K, 


SL 


dL 


m 


section lift coefficient 

flapping damper coefficient, lb-ft-sec/rad (N-m-sec/rad) 

lag damper coefficient, lb-ft-sec/rad (N-m-sec/rad) 

aerodynamic drag moment coefficient in lead-lag equation of motion at given 
blade azimuth angle 

aerodynamic thrust moment coefficient in flapping equation of motion at 
given blade azimuth angle 

incremental drag force on blade element, lb (N) 

offset of center line of flapping hinge from center line of rotor shaft, ft (m) 

offset of center line of lead-lag hinge from center line of flapping hinge, 
ft (m) 


offset of center line of lead- lag hinge from center line of rotor shaft, 
ej + e 2 , ft (m) 


gravitational acceleration, ft/sec2 (m/sec2) 

mass moment of inertia of hub arm e£ about flapping hinge, 


y * m(r - ej^dr, slug-ft^ (kg-m2) 


mass moment of inertia of blade about lead-lag hinge, 



flapping spring constant, lb-ft/rad (N-m/rad) 
lag spring constant, lb-ft/rad (N-m/rad) 
incremental lift force on blade element, lb (N) 
mass of blade per unit of radius, slugs/ft (kg/m) 
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m a 


M d 


Mj 


M S 


M W 


M W 


M W,e 


M (x,^) 

dN 


dP 


r 


r 


c 


r F 


R 

t 

u 
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aerodynamic-force moment, lb-ft (N-m) 
damper moment, lb-ft (N-m) 
inertia-force moment, lb-ft (N-m) 
spring- restraint moment, lb-ft (N-m) 
weight moment, lb-ft (N-m) 

pR 

weight moment of blade at zero angular displacement, j mg(r - e^dr, 
lb-ft (N-m) * 

weight moment of hub arm e 2 at zero angular displacement, 



local Mach number U/a at a given radial station and azimuth angle 

resultant aerodynamic force on blade element perpendicular to blade lag 
plane, lb (N) 

component of aerodynamic force on blade element perpendicular to blade- 
span axis and in the blade lag plane, lb (N) 

distance measured along blade from axis of rotation to blade element, ft (m) 

radius of blade root pocket cutout (i.e., radius at which lifting surface of 
blade begins), ft (m) 

radial distance from flapping-hinge Yp-axis to mass element, measured in 
plane perpendicular to hinge axis, ft (m) 

total radius of undeflected blade, ft (m) 

time, sec 

nondimensional resultant velocity at local blade element, U/ £2R 
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U resultant velocity perpendicular to blade-span axis at blade element, 

ft/sec (m/sec) 


Up component at blade element of resultant velocity perpendicular to both the 

blade-span axis and Up, ft/sec (m/sec) 


Up 


V 

w b 


W e 

x 


component at blade element of resultant velocity perpendicular to the blade- 
span axis and in the Xl-Yl plane, ft/sec (m/sec) 

velocity along flight path, ft/sec (m/sec) 

r R 

blade weight, \ mg dr, lb (N) 

Je t 

P e t 

weight of hub arm e2, \ mg dr, lb (N) 

Je l 

ratio of blade- element radius to rotor-blade radius, r/R 


Xj = e^R 

x 2 = e 2 /R 



X e ,Y e ,Z e 

x f ,y f ,z f 

x l ,y l ,z l 

X S ,Y S ,Z S 


z tip 

“r 


earth reference axes 
flapping axes 
lag axes 
shaft axes 

vertical displacement of blade tip from X s -Y s plane, ft (m) 

blade-element angle of attack, measured from line of zero lift, deg 

angle between shaft axis and plane perpendicular to flight path, positive when 
axis is pointing rearward, deg 
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p 


Kb 


PiP 

^(90-270) 

y' 

6 3 

W b R 2 
6 glh 


blade-flapping-hinge angle at particular azimuth position, rad 
first and second derivatives of /3 with respect to time 
first and second derivatives of j8 with respect to azimuth angle xp 
blade-flapping-hinge angle in azimuth range from 90° to 270°, rad 
mass constant of rotor blade, 

Ih 

blade-flapping-hinge cant angle, deg 


blade leading angle at particular azimuth position, rad 


C,C 


first and second derivatives of £ with respect to time 




77 = 


M^R 

g J h 


first and second derivatives of £ with respect to azimuth angle \f/ 


„ *V R 


© 

e 

9 i 

\ 

A - I e /Ih 

M 

v 


rotor pitch attitude referred to earth axes, deg 

instantaneous blade- section pitch angle, Aq - Aj cos \p - sin \p + 0jx, deg 

difference between root and tip geometric pitch angles, positive when tip 
angle is larger, deg 

V sin <Xq - v 

inflow ratio, 

OR 


V cos o' 

tip- speed ratio, — ~ — 

induced inflow velocity at rotor, positive downward, ft/ sec (m/sec) 
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p mass density of air, slugs/ft® (kg/ m^) 

o rotor solidity, bc/wR 

4> inflow angle at blade element in plane perpendicular to blade- span axis, 

1 • 

tan uf’ deg 

<£ rotor bank angle referred to earth axes, deg 

ip blade azimuth angle, measured from downwind position in direction of rota- 

tion, deg or rad as indicated 

O rotor angular velocity, dip/dt, rad/sec 

Subscripts: 

o initial condition 

F referred to flapping axes 

L referred to lag axes 

max maximum 

s referred to shaft axes 

x,y,z components in x-, y-, and z-directions, respectively 

METHOD OF ANALYSIS 

The equations of motion for a blade with full articulation are derived by equating 
the moments about the hinge points to zero. The general moment expression for the two 
hinge points considered herein is as follows: 


where 

Mj + + M w + Mg + M d = 0 

(1) 

Mj 

inertia moment 


m a 

aerodynamic-force moment 
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I 



weight moment 


M W 


Mg 


spring- restraint moment 


M d 


damper moment 


The two differential equations required are obtained by writing the expressions for 
the individual moments of equations (1). 


Reference Axes 

The basic hub system considered herein is illustrated in figure 1. Included also 
are the axes used to derive the equations of motion. The basic reference axes are the 
shaft axes, which consist of a nonrotating, right-hand, orthogonal system with the 
Z s -axis coincident with the shaft axis and the X s -axis in the plane of the velocity vec- 
tor V, and the Z s -axis. The Y s -axis is perpendicular to the X s - and Z s -axes and is 
positive in the direction of xj/ = 90° on the advancing side of the rotor disk (see fig. 1). 

The additional axes are the flapping axes and lag axes which are rotating, right- 
hand orthogonal axes with origins at the flapping hinge and lagging hinge, respectively. 
The flapping axes consist of an Xp-axis perpendicular to the flapping- hinge axis and in 
the X s -Y s plane, a Yp-axis coincident with the flapping- hinge axis, and a Zp-axis paral- 
lel to the Z s -axis and perpendicular to the Xp- and Yp-axes. The lag axes consist of 
an XL-axis coincident with the hub arm e 2 , a Yp-axis perpendicular to the Xp-axis and 
the lag-hinge axis, and a Z praxis coincident with the lag-hinge axis. It should be noted 
that the lead-lag angle £ is measured in the Xp-Yp plane and that the flapping angle 0 
is the angle between the Xp-axis and the line formed by the intersection of the blade 
plane Xp-Yp and the Xp-Zp plane. 

Lead- Lag Degree of Freedom 

The lead-lag equation- of-motion expressions for moments acting about the lag 
hinge are given in the following sections. 

Inertia moment .- The inertia- moment expression may be written as 

r R 

M I; p = -J m(r - e t )a L dr (2) 

where a L is the component about the lag hinge of the total acceleration of a mass 
element. 

A detailed derivation is not presented herein; however, the procedure used was to 
write the expression for the displacement of a mass element in terms of the shaft- axes 
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coordinates and to obtain the first and second time derivatives of the expression. Then 
by use of appropriate coordinate transformation matrices, the total acceleration may be 
expressed as acceleration components referred to the lag axes. These components may 
then be resolved into the desired resultant component perpendicular to the blade span. 
This procedure yields the following expression for a^: 


t L = ?(r - e t ) - 2^[(r - e t ) cos 2 ^ + 6g)sin 0 + e 2 cos 63 cos(£ + 63) sin /3J 
+ 0 2 |jr - e^sin^ + + 63) + 63 cos 63 sin(£ + 63^ 

+ fi 2 ^e^ + eg cos 0)cos 63 cos 0 sin(£ + 63) - e t sin 63 cos(£ + 63^ 

- fi 2 ^r - e t )sin^ + 63^005^ + 


(3) 


Substituting equation (3) into equation (2) and integrating yields 


M 


. M 


j L = -I h C + 2^2/31^ cos*^ + 63^sin 0 + 2^2/3 63 cos 63 cos(£ + 63^sin 0 


- jS^Ijj sin(£ + 6 3 )cos(£ + 63) - 0 2 e2 cos 63 sin(£ + 63) 


t 

9 Mw r . , 

- £2 “g^ ( e l + e 2 cos P) cos 5 3 cos i 3 sin (? + ^3) “ 

+ sin^ + 63^cos^ + 63) sin 2 0 


e t sin 63 cos^ + 


63)] 


(4) 


Aerodynamic-force moment .- The aerodynamic-force moment about the lag hinge 
is obtained by resolving the element lift and drag forces into a resultant force about the 
lag hinge. It should be noted that the aerodynamic forces are obtained by using a quasi- 
static, independent- element analysis. Figure 2 illustrates the forces acting on the ele- 
ment. The element lift force is 

dL = | pU 2 c z c dr (5) 

and the element drag force is 

dD = | pU 2 c d Q c dr (6) 

The desired resultant force dP is given by 


dP = dD cos <p - dL sin cp 


(7) 
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From figure 3, the aerodynamic-force moment is 



Substituting equations (5), (6), and (7) into equation (8) yields 


( 8 ) 


M 


A,L 



pR 1 9 

e t )sin (p dr - j - pUc^ Q c ( r - e^cos <p dr 


(8a) 


It is assumed that r c is equal to or greater than e^. so that no aerodynamic- moment 
contributions are included for the blade spar inboard of the root pocket cutout. In addi- 
tion, the integration of the lift integral contains the tip- loss factor B so that lift forces 
beyond BR are neglected. 

The additional expressions required to evaluate equation (8a) are as follows: 

U 2 = U T 2 + Up 2 (9) 

where 


U T = S2(r - e t )cos /3 + Oe t cos £ cos )3 + fie^ cos 63 cos(£ + 63)^1 - cos /3) + £(r - e t ) 

+ ORpJcos^ - 63^sin^ + 63)005 0 + sin(i// - 63)003^ + 

+ £2RA sin^ + 63)sin j3 (9a) 


and 


Up = J2RA cos j3 - S2R/J. cos ( '\]s - 63)sin /3 - 0(r - e^sin(? + 63)sin 0 


r - e t )cos(£ + 63) - ^2 cos 63 - S2e t sin 63 sin /3 

(9b) 

-1 U P 
<p - tan 1 

Ut 

(10) 

otp = 9 + (p 

(11) 

Plrt 

11 

2 

(12) 


Expressions (11) and (12) are used to determine the section lift and drag coeffi- 
cients required in equation (8a). The lift and drag coefficients are functions of section 
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angle of attack and Mach number; thus, compressibility and two-dimensional stall effects 
may be taken into account for arbitrary airfoil sections. 

It should be noted that the section angle of attack (eq. (11)) and the velocities Up 
and Up are derived with respect to the lag axes, and consequently, the pitch change 
due to the flapping and lag position is taken into account. In order to decouple the blade 
pitch from the flapping and lag displacement, the following expression may be used for 
the section angle of attack: 

a r = 8 + (p + A6 (13) 

where 

A 9 = tan” * (sin £ tan j3) (13a) 

Weight moment.- The expression for the in-plane weight moments is obtained by 
first orienting the gravity vector with respect to the shaft and then determining the com- 
ponent about the lag hinge from the appropriate coordinate transformations. As illus- 
trated in figure 4, the gravity vector must be transformed through the angles 0 and 
to obtain components referred to the shaft axes. The expression for the transformation 
is 


c ^ 


— 


- 


r -\ 

Sx 


cos © 

0 

-sin 0 


0 

( gy 

> = 

sin $ sin © 

cos $ 

sin $ cos 0 

\ 

; 

0 

g z 


sin © cos $ 

-sin <£ 

cos $ cos 0 

! 

j 

-g 

V- J 

' S 

. 


_ 


J 


By appropriate coordinate transformations, the gravitational acceleration about the lag 
hinge may be shown to be 

gL = g sin(? + 6g)sin 0 cos <i> cos © 

- gjcos(i// - 6g)sin(£ + 6g)cos /3 + sin^ - 5g)cos(£ 

+ gjsin^ - 5g)sin(£ + 6g)cos /3 - cos(^ - 6g)cos^ 

The expression for the weight moment about the lag hinge is 

r R 

M W,L = J et m g L ( r ” e t) dr ( 16 > 


+ Sgjjsin 0 

+ 6gjj sin <3> cos © (15) 
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Substituting equation (15) into equation (16) and integrating yields the following expres- 
sion for the weight moment: 

= + 63^ sin ^ cos $ cos© 

- Jcos^ - f>3)sin(C + 63) cos 0 + sin^tf/ - 63)008 + 63^ sin 0 
+ fsin^ - 6 3 )sin(£ + 6 3 )cos /3 - cos (ip - 6 3 )cos(£ + 6 3 jjsin $ cos ©| (17) 


Spring- restraint moment .- The spring moment which is included in the in-plane 
equation of motion may be used to represent either an actual spring restraint or to 
simulate the structural stiffness of a hingeless rotor. 

The expression for the moment is 

M S,L = _K SL ? ( 18 ) 


where KgL is the lag spring constant. 

Damper moment. - The damper- moment expression has been included to represent 
either a viscous damping moment or simulated structural damping. 


The moment is defined as 




(19) 


where L is the lag damper coefficient. 

Lead- lag equation of motion .- The expressions for the moments about the lag hinge 
have been derived in the preceding sections. Substitution of the five terms (eqs. (4), 

(8a), (17), (18), and (19)) into equation (1) yields the following equation of motion for the 
lag degree of freedom: 


-Ih? + Ijjfi 2 sin(£ + 6 g)cos(? + 6g)sin 2 /3 - I[jj3^sin(iJ + 6 g)cos^ + 63 ) + 21^0/3 cos 2 (? + 63 ) sin p 

n 2 Mw 2 2 

^ cos 63 cos /3 sin(£ + 63) + SI' e t sin 63 cos(£ + 63) — £ 2 ' e 2 cos 63 cos £ sin(£ + 63) 

t t 

j3 2e 2 cos 63 sin(£ + 63)+ 2 fi0e 2 cos 63 cos^ + Sg^sin )3 - M W jjc° s (^ - 6 3 )sin(£ + 6 3 )cos 0 
+ sinfo - 6 3 )cos(^ + 63 jj sin 0 - jsin(^/ - 6 3 )sin(£ + S 3 )cos 0 - cos [if/ - 6 3 )cos(C + 63)j sin $ cos 0 

- sin(£ + 63) sin )3 cos $ cos 0 ^ + J i pU 2 Cj c sin 0 ( r " e t) dr “ § \ pU 2 c d o C cos *(r - e t )dr - Kg L S - = 0 (20) 
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Equation (20) is the dimensional equation of motion. It may be nondimensionalized 

O 

by dividing by -1^ and substituting the following expressions: 


^ = dt,£/ d£ - 

^ dt d xp s 

•• 9 = 


/3 = f2j8 


p = n 2 /3 


e^ = XjR 


e 2 = x 2 R 


e t = x t R 
r = xR 


r c = x c R 


U = uS2R 


P£ R = y . 

h 


W b R 

g J h 


M^R 


= e 


= t] 


J h 


M W.e R _ 


g!h 


= r ?e 


The nondimensional lead-lag equation of motion is 

£ - [sin 2 /3 - / 3 ^jsin(£ + 6g^cos(? + 6 ^- 2 ^ cos^£ + 63^ sin /3 

+ ?7^Xj + Xg cos / 3 )cos 6g cos /3 sin(£ + 63^ - rjx^ sin 6g cos^ + 6g) 
-2 

+ 77/3 x 2 cos 63 sin(£ + 63) - 277/8x2 cos 6g cos(£ + 6g)sin 0 

+ -^|cos(* - 6^sin^£ + 63^ cos 0 + sin^ - 6gjcos(£ + Sgjjsin 0 

- jsin^ - Sgjsin^ + Sg^cos 0 - cos(x{/ - Sg^cos^ + 6gJJsin 3 > cos 0 

(C + ^3) sin p cos $ cos 0 + C Md + ? = 0 

n Ih^ 


- sin 


where 


P 1*0 1 9 pB i Q 

C M D = 1 2 r ’ U C d,o cos H x - x t) dx “1; 2 r ’ U C * Sin 0 ( x • *$** 

x c 


( 21 ) 


( 22 ) 


(22a) 


Flapping Degree of Freedom 

The equation of motion for the flapping degree of freedom is obtained by again 
determining the expanded moment expressions of equation (1). The expressions for 
moments about the flapping hinge are given in the following sections. 
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Inertia moment.- The rigid- blade inertia moment about the flapping hinge may be 


written as 




R 


mr F a F dr 


(23) 


where r p is the radial distance from flapping- hinge Yp-axis to mass element, mea- 
sured in the plane perpendicular to the hinge axis. 

By using the same procedure as was described for the acceleration about the lag 
hinge, the acceleration a F and the moment arm r F may be expressed as 


a F = 0(r - e^cos 63 + £2 2 ^r - e^sin /3 cos p + S2 2 ej cos 63 sin (3 ^ e^ 


r F = ( r “ e l) cos 6 3 
and 


( e l = r = e t) 


(24a) 

(24b) 


a F = n2 ( r " e t) c0S (^ + 6 3) sin P cos P + £( r " e t) cos (^ + 6 3) “ 2 ^( r " e t) sin (^ + 6 3) 

• •• 2 

+ 2 S 2 i^r - e^cos^ + Sg^sin /3 + /Seg cos 63 + O e^ cos 63 sin (3 

2 

+ S 2 e 2 cos 6g sin j8 cos /3 
r p = e 2 cos 63 + (r - e^cos^ + 63) 


! t = r = R ) 

(25a) 

IT 

VII 

u 

VII 

-M 

(25b) 


By substituting equations (24) and (25) into equation (23) and integrating with appro- 
priate limits, the inertia moment may be expressed as 

M if = -I h n 2 cos 2 (£ + 6 3 )sin j3 cos j3 - I h /3 cos 2 (C + 63 ) 


2I h /3? sinp + 6 3 Jcosp + 63 )- 2l h S2£ cos 2 JC + Sg)sin /3 


M W 


M 


W 


- 2-^- egfi cos 63 cos^ + 63 ^ sin j3 cos /3 - 2 — — e 2 /3 cos 63 cos^ + 63 j 


M 


W 


M W 


ejO cos 63 cos^ + Sg^sin j3 + 2 — - e 2 /3£ cos 63 sin^ + 6 gj 
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Mttt . W, 9 2 2 

- 2 —^r- e 2 &£ cos 6g cos^£ + 6^ sin /3 - -^- e 2 S2 cos 6g sin ]3 cos /3 

W. n 2 W b 2 2 

- -g-2 e 2 0 cos 6g - -g^ ejegO cos 6 3 sin ^ 

f 

- I e fi 2 cos 2 6 3 sin /3 cos )3 - I e /3 cos 2 6g - — ^ e 1 fi 2 cos 2 6g sin £ (26) 


Aerodynamic-force moment .- The expressions for the aerodynamic forces on a 
blade element are given in equations (5) and (6). From figure 2, the component of these 
forces (i.e., lift and drag) normal to the blade lag plane is 

dN = dL cos <p + dD sin (p (27) 

The aerodynamic-force moment is therefore 

M A,F = I r R r F dN < 28 > 

r c 

where 


r F = e 2 cos 6g + (r - e^cos^C + 6g) 

By substituting equations (5), (6), and (27) into equation (28), the aerodynamic-force 
moment may be expressed as 

1 A F = J \ c cos 4>je 2 cos 6g + ^r - e^cos^ + 6gj[Jdr 


(28a) 


M, 


§ ^ pU 2 c d o c sin ^[ e 2 COS 5 3 + ( r “ e t) cos (^ + 6 3)J dr 


(29) 


Weight moment. - The components of the gravitational- acceleration vector referred 
to the shaft axes are given in equation (14). By appropriate coordinate transformations, 
the gravitational acceleration about the flapping hinge may be shown to be 


gp = -gj^sin 0 cos^p - 6g^ - sin <3> cos 0 sin^ip - 6gjj sin /3 - g cos $ cos 0 cos /3 (30) 

The weight moment about the flapping- hinge axis, defined as positive in the direc- 
tion of positive j3, is 
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f R 

M W,F = 1 m ^F r F dr 
e l 


( 31 ) 


where 


F = ( r - e i) cos 6 3 

'lb' 

IIA 

IIA 

F = e 2 cos 6 3 + ( r ' e t) C0S P + 6 3) 

'pT 

VII 

u 

VII 

-M 


> (31a) 


By substituting equation (30) into equation (31) and integrating with appropriate limits, 
the weight moment may be expressed as 

M W F = -(< cos^C + 63) + e + W^eg^cos 6^j jcos $ cos© cos ft 
+ sin 0 sin 0 cos(\f/ - 6^ - sin <f> cos © sin /3 sin^ - 6gjj 


(32) 


Spring- restraint moment .- The spring moment about the flapping hinge is included, 
as was done for the lag hinge, to permit simulation of structural stiffness or an actual 
spring at the hinge point. The moment expression is simply 


M 


S,F ~ ^SF^ 3 


(33) 


Damper moment. - The damper moment is included also for generalization of the 
rotor system. The moment expression is 


M d,F = ” C d,F^ 


(34) 


Flapping equation of m otion.- By substituting the five moment expressions derived 
in the preceding sections (i.e., eqs. (26), (29), (32), (33), and (34)) into equation (1) and by 
using the nondimensionalizing parameters given in equations (21), the flapping equation 
of motion may be shown to be 

jsjcos 2 (? + + eXg^cos^Sg + 27yxg cos 6 3 cos^ + 6 3 jj + |cos 2 (? + 5 ^ 

+ + exg^jcos^ + 277Xg cos 6g cos^? + 6^ sin j8 cos /3 

+ IjeXjXg + 77 e Xj^cos 2 6g + r/x j cos 6g cos(? + 6^ + 277 X 2 ? cos 6g cos(? + 


16 



( 35 ) 


+ 2? cos 

. e 


+ 6 3 )J sin P - 2^Jsin^ + 6g)cos(£ + 6g)+ rpt 2 cos 6g sin(? + 6gjj - C Mt 
|r7 cos^ + 6g) + ^n e + ex 2 )cos 6^J jcos $ cos 0 cos 0 + sin 0 sin /3 cos(i// - 6g) 


- sin 4? cos 0 sin /3 sin^ - 6gjl + j3 + — /3 = 0 

J h £2 


where 


Cm-j- = § ^ y’u^c^ cos cos ®3 + ( x ” x t) cos (^ + <ix 


§ ^ y'u 2 c d o sin ^[*2 cos 6 3 + ( x " x^cos^ + 63 |j dx (35a) 


Solution of the Equations of Motion 

The equations of motion (eqs. (22) and (35)) are nonhomogeneous and nonlinear with 
varying coefficients that contain complex integral expressions when such items as stall 
and compressibility are considered. Thus, it is not feasible to obtain an explicit solu- 
tion and a numerical iterative method is required. 

The procedure used to obtain the results presented is outlined as follows: 

1. At an initial azimuth position \p Q , assume initial conditions for £, £, (3, 

and js. 

2. Compute the thrust moment and drag moment at the selected 

azimuth position by evaluating the integrals at specified radial stations and integrating 
the results by numerical methods. 

3. Solve equation (22) for £ and equation (35) for j§. 

4. Repeat steps 2 and 3 at the next azimuth position, using values of £, £, /3, 

and j3 obtained by integrating the previously determined quantities £ and /3. Several 
techniques are available for performing the integration. The method used herein is the 
Runge-Kutta method described in reference 5. 

5. Continue the process at successive azimuth positions until (3 and £ reach 
steady-state values (i.e., repeat for successive rotor revolutions) or one or both of the 
angles diverge. 

The method of solution outlined above yields the transient blade motion from some 
arbitrary initial condition to a steady- state flight condition. The method may also be 
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used to compute the transient response of the blade motion due to control inputs or angle- 
of-attack changes by changing the appropriate input information after the initial steady- 
state flight condition has been achieved. Rotor-performance characteristics and steady 
blade motions may also be obtained from the forces and motions computed when the blade 
motion has reached a steady- state condition. The data presented herein do not cover 
rotor-performance aspects; however, as an aid in checking the equations and programing, 
the performance results obtained with this program have been compared with other com- 
puter programs (e.g., ref. 2) and very good correlation was found. 

APPLICATION AND DISCUSSION 

The applications and results presented herein are intended to illustrate the prob- 
lems associated with rotor operation at high tip- speed ratios and to define limiting con- 
ditions of rotor operation. It must be emphasized that the boundaries presented are for 
a particular rotor-hub configuration and are not intended to be applicable to all design 
configurations (e.g., configurations with added hinge restraint and/or dampers). The 
boundaries presented do, however, indicate the approximate range of tip- speed ratios 
wherein the destabilizing aerodynamic spring forces will have an adverse effect on blade 
motion. 

The following table lists the primary geometric parameters of the rotor hub and 
blade which were used in the examples presented herein: 

a 0.1 

0 

0.05 

0.2 

0 

0 

0 

0 

0 

0 

0.97 

Varied as indicated 

Constant (value compatible with y') 

It should be noted that the in-plane damper has been omitted, and consequently, the 
lag motion will exhibit only the small amount of aerodynamic damping. 

The airfoil section characteristics used herein are those given in reference 6. 

The section characteristics include compressibility effects; however, in order to 


X 1 ' 
x 2 . 

x c . 

C d,F 

K SF 

C d,L 

K SL 

63 . 

*1 * 
B . . 

/. • 
m . . 
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minimize such effects on the results presented herein, the rotor tip speed was adjusted 
to maintain an advancing tip Mach number of 0.8. 

Significance of the Angles (3 and £ 

The hub restraint specified in the derivation of the equations of motion (i.e., lag 
hinge at or outboard of the flapping hinge) requires the use of a flapping angle f3 which 
does not directly describe the physical location of the blade. It does, however, give the 
angular displacement of the flapping- hinge bearing. For example, at extreme lag 
angles 5, it is possible to attain a 90° flapping angle with the blade only slightly out of 
the hub plane. This may be more clearly illustrated by writing the expression for the 
blade-tip displacement with respect to the hub plane. From figure 1 the expression for 
coincident hinges at the center of rotation is simply 

z tip = ^ cos ^ s * n ^ 

It is apparent that as £ approaches 90° (cos £ — 0), the tip displacement 
approaches zero regardless of the amplitude of /3. Thus, at extreme lag amplitudes, 
the motion of the flapping hinge must be large for even moderate tip displacement. Since 
this unique condition exists only at large amplitudes of lag, the flapping angle can gen- 
erally be considered as a measure of blade motion out of the hub plane for small lag 
angles (e.g., when cos £ « 1). 


Stability Examples 

Some preliminary applications and results obtained from a numerical solution of 
the equations of motion derived herein are presented in reference 7. Figure 5 herein, 
taken from reference 7, shows the blade-motion stability boundaries developed by intro- 
ducing the same initial flapping- angle displacement at two different initial azimuth posi- 
tions and solving for the transient blade motions of an unloaded (nonlifting) rotor. The 
two boundaries illustrate the sensitivity of blade- motion stability to the azimuth position 
at which a displacement disturbance occurs. It should be noted that the boundaries pre- 
sented are essentially absolute limits, in that they represent conditions where hinge- 
motion amplitudes exceed 90°. Consequently, acceptable rotor operation (i.e., hinge- 
motion amplitudes less than blade mechanical stops) would be limited to tip- speed ratios 
less than those indicated by the boundaries. In addition, the disturbance amplitude intro- 
duced is relatively large and, thus, the lower boundary is in all likelihood conservative. 
These boundaries, however, were developed for a nonlifting rotor system and it would be 
expected that the addition of a lifting load would aggravate the condition because lift 
produces a destabilizing moment. An example of the shift which can occur for a lifting 
rotor is presented in reference 7. 
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The fact that stability boundaries may be established over a range of tip- speed 
ratios for various amplitudes and phasing of the disturbance is significant because it is 
possible to operate a rotor beyond a specific boundary and then encounter a disturbance 
of sufficient magnitude to produce a divergent condition. 

As a measure of the sensitivity of the boundaries to the initial flapping- angle input, 
the variation of the limit established for a mass constant y' of 1.6 for varying initial 
flapping angles is presented in figure 6. The boundary presented was established in the 
same manner as the boundaries in figure 5. That is, the blade was released from the 
initial flapping- angle amplitude 0 O as indicated and the transient solution of the equa- 
tions of motion was obtained. As indicated in figure 6, the increase in the tip- speed 
ratio at which unstable motion occurs for decreasing /3 0 is relatively linear for the 
range ]S 0 = 0.10 to /3 0 = 0.20. Below /3 0 = 0.10, the boundary moves more rapidly to 
higher tip- speed ratios, as might be expected for such small disturbances. Again it 
should be emphasized that these boundaries are, in effect, absolute limits since the blade 
amplitudes which occur in the stable region near the boundary are unacceptable in terms 
of realistic structural limits. 


Limited- Response Boundary 

An example of the amplitudes obtained for stable transient motion of the blade near 
the boundary for /3 0 = 0.20 (presented in fig. 6) is shown in figure 7. The transient 
motion presented is characteristic of the response obtained when the blade is released 
with an initial displacement /3 0 at \p Q - 90°. That is, the flapping angle increases as 
a result of the negative effective aerodynamic spring in the forward quadrants and the 
blade swings forward (i.e., leads) as a result of the Coriolis forces resulting from the 
positive flapping velocity. 

The response curve also illustrates another trend which is generally characteristic 
of the motion for the prescribed initial conditions: for stable conditions near the bound- 
ary, the maximum flapping amplitude occurs within the first half revolution after the blade 
is released. As was pointed out in reference 7, this increase in flapping amplitude occurs 
because the negative (destabilizing) aerodynamic spring forces in the forward quadrants 
are larger than the restoring centrifugal spring force. Thus, a more realistic boundary 
in terms of acceptable rotor operation could apparently be based on the tip-speed ratio at 
which the restoring spring forces are at least equal to the negative spring forces during 
the first half cycle after the blade is released at 90° azimuth. Thus, upon release of the 
blade from an initial amplitude, the flapping angles do not exceed the initial amplitude 
(i.e., j8 s p o ) for the first half cycle. A boundary based on this criterion is presented 
in figure 8, and represents the approximate tip- speed ratio at which the rotor blade 
encounters regions of neutral stability in the forward quadrants. Also included is the 
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absolute limit of figure 5. The area between the two boundaries may be considered as 
an area of marginal stability. That is, it is an area in which the sensitivity to distur- 
bances is increasing rapidly as a result of the contribution of the destabilizing aero- 
dynamic spring forces in the forward quadrants of the rotor disk. 

An example of the blade transient motion which is characteristic of conditions 
along the lower boundary of figure 8 is presented in figure 9. The significant trend of 
the transient motion is the relatively flat slope of the flapping amplitude for the first 
part of the initial revolution as opposed to the large increase shown for the condition 
presented in figure 7. This relatively flat slope indicates that the nonlinear centrifugal 
and aerodynamic spring forces are approximately equal for the conditions imposed. In 
addition, the remaining revolutions indicate that the flapping motion is rapidly damped 
and that, as expected, the in-plane (lagging) motion has very little damping. The lagging 
amplitude, however, is not large because the coupling terms (Coriolis forces) are small 
for the conditions presented. 

Since the boundary presented in figure 8 is based on somewhat arbitrary initial 
conditions, it is of interest to illustrate the significance of the boundary. This was done 
by establishing a boundary based on a criterion that the maximum flapping does not 
exceed 15° when a nonlifting rotor is subjected to a 30-ft/sec (9.14-m/sec) vertical 
gust. The vertical gust was introduced as a step function in the inflow ratio X. The 
resultant boundary is presented in figure 10 and compared with the boundary of figure 8. 

It should be noted that the gust input was made with the blade at an initial azimuth angle 
of zero degrees with zero initial flapping displacement. The agreement between the two 
boundaries is surprisingly good in view of the fact that they are based on two diverse 
criteria. In addition, it appears that the boundary based on the criterion of neutral 
spring forces in the forward quadrants ( ^(90-270 ) “ ^°) £* ves a ver Y good indication of the 
tip -speed ratio at which blade- motion stability begins to deteriorate. Thus, the sensi- 
tivity to disturbances would probably be of concern when operating at tip- speed ratios 
above the indicated boundary. 


CONCLUDING REMARKS 

The equations of motion for a fully articulated rotor system (flapping and lead- lag 
degrees of freedom) have been derived and a method for obtaining numerical solutions 
of the equations is presented. The equations are applicable to a wide range of rotor 
variables (e.g., airfoil section, twist, root cutout, and hinge geometry) and can account 
for two-dimensional stall and compressibility effects. The hinge restraints considered 
require that the lag hinge be coincident with or outboard of the flapping hinge. With this 
hinge restraint, the flapping angle j8 describes the flapping- hinge motion and does not 
indicate directly the blade motion with respect to the hub plane. 


The equations were used to define blade motion stability boundaries in terms of an 
absolute limit (i.e., amplitudes greater than 90°) and to define a boundary based on the 
criterion of neutral spring forces in the forward quadrants when the blade is released 
with an initial flapping- angle displacement at an initial azimuth angle of 90°. The latter 
boundary is thus a measure of the tip- speed ratio at which the destabilizing aerodynamic 
spring forces in the forward quadrants are approximately equal to the stabilizing centrif- 
ugal forces. Consequently, rotor operation above the boundary could be considered to be 
in an area of marginal stability in the sense that the rotor sensitivity to disturbances is 
increasing rapidly. A comparison of the boundary based on the neutral- spring criterion 
with a boundary based on a maximum flapping- amplitude limit of 15° for a 30-ft/sec 
(9.14 m/sec) vertical gust indicated good correlation. 
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Figure 7.- Transient motion near unstable boundary. V/QR = 1.4; y‘ = 1.6. 



Figure 8.- Comparison of boundary based on neutral spring forces in forward quadrants with absolute- limit boundary. 

= 7t/2 rad; J3 0 = 0.2 rad. 
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Figure 10.- Comparison of neutral-spring-force boundary with blade-flapping-amplitude boundary. 
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